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Summary 

Under  USAF  contract,  Microsystems  and  Nanotechnologies  performed  research  and  technology 
development  in  three-dimensional  (3D)  nanobiocomputing  architectures.  The  major  accomplishments  are: 

1.  A  novel  design  concept  for  Integrated  Circuits  (ICs)  and  computing  architectures  in  3D  space  was 
developed  and  demonstrated.  The  proposed  design  is  supported  by  novel  software  and  envisioned 
hardware  technologies.  For  complex  ICs,  we  performed  the  design  by  utilizing  neuronal 
hypercells  (N  -hypercells)  within  a  neuronal  aggregated  network.  At  the  system  level, 
benchmarking  ICs  (that  have  been  widely  used  to  evaluate  and  assess  2D  Very  Large  Scale 
Integration  /  Ultra  Large  Scale  Integration,  VLSI/ULSI,  design)  were  examined  as  proof-of- 
concept  3D  ICs  for  the  envisioned  biomolecular  computing  architectures. 

2.  Neuronal  hypercells  that  perform  complex  logics  were  designed.  The  aggregated  N -hypercells 
form  neuronal  3D  N -hypertopologies,  thereby  implementing  nanobiocomputing  architectures. 
Interconnected  functional  multi-terminal  biomolecules  form  molecular  gates  (Mgates)  that 
comprise  3D  N -hypercells. 

3.  Representative  Compute  Aided  Design  (CAD)  components  were  developed  to  coherently  support  the 
reported  design  tasks  for  3D  ICs. 

4.  A  DNA  transistor  (DNAT),  as  a  possible  proof-of-concept  three-terminal  hybrid  bioelectronic  device, 
was  modeled,  simulated  and  analyzed. 

We  departed  from  the  two-dimensional  (2D)  paradigm  and  2D  topology/organization/architecture 
of  ICs  by  developing  a  novel  technology  that  utilizes  a  three-fold  innovative  solution: 

•  hardware  (novel  3D  molecular  devices,  N -hypercells,  modules  and  systems); 

•  software  (technology-specific  CAD  supported  by  new  synthesis  and  design  methods); 

•  nanotechnology. 

Thus,  in  general,  we  propose  a  revolutionary  3D1 2 3  (Hardware-Software-Nanotechnol°gy)  technology  that 
ensures  super-high-performance  computing  and  information  processing.  We  focused  on  development  of  a 
feasible,  practical,  affordable  and  superior  technology  for  massive  parallel  distributed  computing  and 
information  processing.  The  proposed  technology  will  allow  one  to  perform  design  coherently  integrating 
synthesis  tasks  at  the  device/gate  and  system  levels.  We  depart  from  conventional  2D  VLSI  and 
Complimentary  Metal  Oxide  Semiconductor  (CMOS)  paradigms.  A  novel  3D3  technology  will  offer: 

1  Enormous  military  advantages  because  newly  designed  super-high-performance  platforms  will 
guarantee  information  processing  preeminence,  computing  superiority  and  memory  supremacy; 

2  Very  strong  commercial  potential  with  immediate  applications  in  design  of  new  generations  of 
preeminent  processors  and  memories; 

3  Sound  technology  transfer  feasibility  to  future  Air  Force  systems. 
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1.  Introduction 


This  contract  started  December  27,  2004.  The  overall  objective  was  to  devise,  design  and 
demonstrate  novel  super-high-performance  three-dimensional  (3D)  nanobiocomputing  architectures  with 
aggregated  3D  neuronal-hypercells  (K  -hypercells)  that  integrate  interconnected  functional  multi-terminal 
electronic  nanobiodevices  (biomolecules).  The  specific  objectives  were  to: 

1  Devise  and  demonstrate  3D  nanobiocomputing  architectures  designed  utilizing  a  network  of 
aggregated  X  -hypercells  (modeling,  analysis  and  software  development)  that  form  ICs. 

2  Design  X -hypercells  that  perform  complex  logics.  Develop  a  library  of  basic  X -hypercells. 

3  Develop  representative  components  of  3D  nanobiocomputing  CAD  with  a  focus  on  X  -hypercells. 

4  Design  8-bit  arithmetic  logic  unit  (ALU)  as  a  proof-of-concept  platform  to  prove  and  demonstrate 
super-high-performance  computing  capabilities  of  3D  computing  architectures. 

5  Model  and  analyze  DNA-based  transistors  (DNAT).  Examine  the  electronic  properties  of 
biomolecules  as  electronic  nanobiodevices  (modeling,  simulation  and  software  developments). 

6  Analyze  the  system-  and  device-level  integration,  aggregation,  compliance  and  compatibility. 

The  research  at  the  system  and  device  level  was  focused  on  the  solution  of  a  number  of  major 
problems.  The  technology  developments  were  concentrated  on  the  following  four  major  tasks: 

1.  Devise  and  demonstrate  a  novel  3D  super-high-performance  nanobiocomputing  architecture 
utilizing  X -hypercells.  Developing  the  unified  bottom-up  and  top-down  design  taxonomy, 
aggregate  3D  X  -hypercells  within  3D  X  -hypertopologies.  Design,  demonstrate,  analyze  and 
evaluate  a  proof-of-concept  ALU. 

2.  Design  and  Integrate  X -hypercells  that  perform  complex  logic.  Develop  a  library  of  basic  X- 
hypercells  for  3D  computing  architectures.  The  major  problems  to  be  addressed  are:  synthesis  and 
logic  design,  as  well  as  X  -hypercell  aggregation. 

3.  Develop  representative  components  of  3D  nanobiocomputing  CAD  with  a  focus  on  3D  ICs 
comprised  from  X -hypercells  that  integrate  functional  multi-terminal  electronic  nanobiodevices. 
Solve  the  following  problems:  (1)  logic  design  and  analysis;  (2)  topology  synthesis  and 
aggregation;  (3)  steady-state  and  dynamic  analysis  of  DNAT. 

4.  Develop  high-fidelity  and  compact  (tractable)  mathematical  models  for  DNAT  examining 
electronic  properties.  These  models  must  support  heterogeneous  simulations  and  data-intensive 
analysis  of  electronic  nanobiodevices  aggregated  in  X -hypercells.  Verify  models  using  the 
experimental  data  for  the  electron  transport  in  DNAT. 

For  the  aforementioned  major  tasks,  in  order  to  accomplish  the  overall  and  specific  objectives. 
Microsystems  and  Nanotechnologies  successfully  accomplished  the  following: 

For  Task  1.  We  devised  and  demonstrated  novel  3D  nanobiocomputing  architectures  utilizing  X- 
hypercells  by  designing  3D  ICs.  The  synthesis  is  performed  using  the  unified  bottom-up  and  top-down 
design  taxonomy  aggregating  3D  X -hypertopologies.  This  task  is  supported  with  programs  developed  in 
the  C  and  MATLAB  environments.  Demonstration,  analysis  and  evaluation  are  performed  for  a  proof-of- 
concept  ALU  and  other  ICs. 

For  Task  2.  We  designed  X-hypercells  that  perform  complex  logic.  We  developed  a  library  of 
basic  X -hypercells  within  nanobiocomputing  architectures  and  3D  X -hypertopologies. 
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For  Task  3.  We  developed  representative  components  of  3D  nanobiocomputing  CAD 
(representative  CAD  tools)  with  a  focus  on  design  of  3D  ICs  comprised  from  X -hypercells.  We 
coherently  supported  the  basic  and  applied  research  developing  MATLAB-  and  C-based  software. 

For  Task  4.  We  modeled,  simulated,  analyzed  and  assessed  DNAT  as  an  electronic 
nanobiodevice  for  X -hypercells. 


2.  Methods,  Assumptions,  and  Procedures 

Proposed  Solution  (Methods)  -  Our  proposed  solution  and  technology  are  based  on  the  design 
of  3D  computing  architectures  and  ICs  utilizing  a  three-fold  solution  (method): 

1 .  innovative  hardware  (device/module/system  levels  and  3D  ICs), 

2.  novel  software  (CAD  supported  by  new  synthesis  and  design  methods), 

3.  nanotechnology-centered  hardware-software  co-design. 

In  particular,  Microsystems  and  Nanotechnologies  proposed  and  demonstrated  a  revolutionary 
3D3  (Hardware-Software-Nanotechnology)  technology.  Our  solution  ensures  information  processing  preeminence  and 

computing  superiority.  By  introducing  the  proposed  3D3  computing  technology,  we  depart  from 
conventional  planar  ICs  topologies  as  well  as  IC  design  techniques,  such  as  VLSI,  ULSI  and  post  ULSI. 

A  technology  of  massively  parallel  distributed  (pipelined  due  to  systolic  processing)  computing 
and  information  processing  in  3D  is  utilized  using  devised  X  -hypercells.  In  general,  X -hypercells 
significantly  increase  the  number  of  bits  of  information  processed  and  exchanged  to  compute  complex 
switching  functions.  The  redundancy  and  parallelism  ensures  robust  fault-tolerant  computing.  The  logic 
design  and  mappings  are  simplified  due  to  the  linear  decision  diagram  and  linear  systolic  array  for 
information  processing  and  manipulation.  We  integrated  fundamental  and  applied  research  within  innovative 
computing  and  molecular  electronics  technologies. 
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3.  Results  and  Discussions 


3. 1.  DESIGN  TAXONOMY  AND  LOGIC  DESIGN  OF  THREE-DIMENSIONAL  ICS 


3.  1 .  1 .  Synthesis  Taxonomy 

To  design  novel  nanobiocomputing  architectures,  we  utilized  an  innovative  unified  top-down  and 
bottom-up  synthesis  taxonomy.  The  proposed  taxonomy  utilized  an  x-design  flow  map  as  reported  in 
Figure  l.a.  At  the  system  and  device  levels,  the  core  design  themes  were  integrated  within  the  following 
four  domains: 


•  Devising  (synthesis)  with  validation;  •  Modeling  -  analysis  -  evaluation; 

•  Design  -  optimization;  •  Nanotechnology  -  fabrication. 


A  top-down  synthesis  was  performed  in  the  proposed  3D  super-large-scale  integration  (SLSI) 
design  in  order  to  synthesize  3D  nanobiocomputing  platforms  implemented  by  aggregating  N -hypercells 
within  ICs.  These  N -hypercells  were  formed  from  molecular  gates  (Mgates)  as  illustrated  in  Figure  Lb. 
The  bottom-up  synthesis  was  focused  on  biomolecular  electronic  devices  aggregated  as  functional 
molecules  that  comprise  Mgates  forming  N -hypercells.  A  unified  top-down  and  bottom-up  synthesis 
taxonomy  was  utilized. 


Devising  a k  Nanotechnology 
Synthesis  ~  /  Fabrication 


(a) 


Multi-terminal  gates 
(electronic  nanodevices 
engineered  from 
biomolecules 


(b) 


Figure  1.  a)  X-design  flow  map; 

b)  K -hypercells  Dijk  formed  utilizing  multi-terminal  Mgates  comprised  from  molecular  complexes 


Figure  2  reports  the  major  tasks  within  the  top-down  and  bottom-up  syntheses  that  were 
concurrently  supported  by  the  envisioned  CAD  to  perform  SLSI. 
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Figure  2.  Nanobiocomputing  architectures  and  associated  design  tasks 


The  neuroarchitectronics  paradigm  serves  as  the  basis  for  the  design  of  novel  super-high- 
performance  biomolecular  computing  platforms  which  can  be  implemented  by  utilizing  biomolecular  3D 
ICs.  Synthesis,  analysis  and  design  of  combinational  biocomputing  hardware  are  performed  with  the 
design  of  3D  ICs.  Using  a  conventional  neuroscience  viewpoint,  we  prototype  a  neuron  as  a  biomolecular 
electronic  module,  while  networked  neuronal  aggregates  are  prototyped  by  X -hypercell  aggregates.  The 
X  -hypercell  aggregates  perform  processing  and  computing.  Figure  3  documents  the  proposed  3D 
nanobiocomputing  architecture.  The  problem  is  to  concurrently  design  3D  ICs. 


Figure  3.  Nanobiocomputing  architecture 
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3.  1.  2 .  Design  of  3D  ICs:  Data  Structure ,  Decision  Diagram ,  and  Hypercells 

An  innovative  solution  to  perform  the  system-level  logic  design  for  3D  ICs  has  been  found  by  the 
Microsystems  and  Nanotechnologies  team  as  reported  in  this  section.  We  departed  from  2D  logic  design 
(VLSI,  ULSI  and  postULSI)  as  well  as  IC  topologies  reported  in  [1].  In  general,  the  proposed  SLSI 
technology  mimics  highly  hierarchical  3D  neuronal  bio  processing  platforms  (brain).  Our  fundamental 
and  technological  solutions  prototype  the  system  architecture,  topology  and  organization  observed  in 
nature.  This  sound  solution  reflects  the  envisioned  device-level  outlook  and  fabrication  technologies. 
Sections  3.2  and  3.4  introduce  and  examine  possible  molecular  electronic  devices  to  implement  N- 
hypercells.  Three  dimensional  architectures/organizations/topologies  and  3D  molecular  electronics, 
unified  by  the  developed  3D  design  concept,  lead  to  a  need  for  novel  CAD  developments.  The  CAD 
solution,  as  reported  in  section  3.3,  should  be  reusable  and  technology-dependent.  The  use  of  3D  N- 
hypercells,  covered  in  section  3.2,  as  baseline  primitives  to  design  computing  platforms,  results  in  a 
technology-relevant  solution.  Our  solution  offers  a  great  deal  of  adaptability  from  the  design  and 
fabrication  technology  viewpoints  ensuring  a  coherent  correspondence  with  bioelectronics  at  the  device 
and  system  levels 

For  2D  CMOS  ICs,  the  decision  diagram  (unique  canonical  structure)  is  derived  as  a  reduced 
decision  tree  by  using  topological  operators.  In  contrast,  for  3D  ICs,  a  new  class  of  decision  diagrams  and 
synthesis  methods  must  be  developed.  We  developed  a  novel  linear  decision  diagram  that  is  mapped  by 
3D  X  -hypercells.  The  hypercell  (cube,  pyramid,  hexagonal  or  other  3D  topological  clusters)  is  a  unique 
canonical  structure  which  provides  a  reduced  decision  tree.  These  hypercells  are  synthesized  by  using 
topological  operators  (deleting  and  splitting  nodes).  We  analyze  optimal  and  suboptimal  (technology- 
specific)  topology  mappings  of  complex  switching  functions.  The  major  optimization  criteria  are:  (1) 
minimization  of  decision  diagram  nodes  and  circuit  terminals;  (2)  simplification  of  topological  structures 
(proposed  linear  arithmetic,  versus  nonlinear  arithmetic,  leads  to  simple  synthesis  and  straightforward 
embedding  of  linear  decision  diagrams  into  3D  topologies);  (3)  minimization  of  path  length  in  decision 
diagrams;  (4)  simplification  of  routing  and  verification  tasks.  The  reported  advantages  result  in  power 
dissipation  reduction,  evaluation  simplicity,  testability  enhancement,  as  well  as  other  important  features. 
For  example,  the  switching  power  is  not  only  a  function  of  devices/gates/switches  (bipolar  junction 
transistors  and  field-effect  transistors,  BJTs  and  FETs,  for  CMOS),  but  also  a  function  of  circuit  topology, 
design  methods,  routing,  dynamics,  switching  activities  and  other  factors  that  can  be  optimized.  We  have 
developed  a  novel  CAD-supported  design  technology  that  allows  one  to  design  high-performance 
nanobiocomputing  architectures.  The  proposed  methods  and  technology,  concurrently  supported  by 
software,  allow  one  to: 

1.  Synthesize  and  aggregate  X  -hypercells  thereby  designing  3D  nanobiocomputing  architectures; 

2.  Perform  logic  design  utilizing  novel  representations  of  data  structures; 

3.  Devise  and  verify  a  new  binary  decision  diagram  synthesis  concept; 

4.  Initiate  the  developments  of  representative  CAD  tools  to  concurrently  support  logic  design. 

The  proposed  concept  will  result  in  SLSI  that  utilizes  a  coherent  top-down/bottom-up  synthesis 
taxonomy  and  nanoarchitectronics.  Novel  3D  nanobioarchitecture  hardware  is  envisioned  to  be 
implemented  as  aggregated  N -hypercells  as  the  fabrication  technologies  (molecular  electronics)  reach 
technical  and  technological  maturity,  economical  feasibility  and  technology  acceptance.  A  number  of 
promising  innovative  paradigms  have  been  explored  in  conventional  microelectronics  [1]  and  emerging 
nanoelectronics  [2-7].  Addressing  the  need  for  fundamental  development  at  the  system  and  device  levels, 
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we  examine:  (1)  design  complexity  (current  CAD-supported  methods  do  not  allow  one  to  design  ICs  with 
a  number  of  gates  more  than  1,000,000);  (2)  device  functionality  (quantum  phenomena  must  be  utilized  in 
molecular  complexes  which  are  fabricated  using  nanotechnology). 

The  binary  decision  diagrams  (BDD)  for  representing  Boolean  functions  are  used  [8].  Reduced- 
order  and  optimized  BDDs  ensure  large-scale  data  manipulations  and  are  used  to  perform  the  logic  design 
and  circuitry  mapping  utilizing  hardware  description  languages  [9,  10].  The  design  scheme  is 

Function  (Circuit)  <->  BDD  Model  <->  Realization. 

The  dimension  of  a  decision  diagram  (number  of  nodes)  is  a  function  of  the  number  of  variables 
and  the  variables  ordering.  In  general,  the  design  complexity  is  0(n3).  This  complexity  significantly  limits 
the  design  capabilities  in  terms  of  ICs  complexity  (number  of  gates)  for  which  the  design  can  be 
performed.  The  commonly  utilized  word-level  decision  diagrams  further  increase  the  complexity  due  to 
processing  of  data  in  word-level  format.  Therefore,  novel  concepts  are  needed.  We  devised  an  innovative, 
sound,  software-supported  design  approach.  In  particular,  novel  methods  in  data  structure  representation 
and  data  structure  manipulation  were  developed  by  the  Microsystems  and  Nanotechnologies  team  to 
ensure  design  specifications  and  objectives.  We  synthesize  and  utilize  the  linear  word-level  decision 
diagrams  (LWDDs)  that  allow  us  to  perform  the  compact  representation  of  logic  circuits  utilizing  linear 
arithmetical  polynomials  (LP)  [11,  12].  The  design  complexity  becomes  0(n).  The  proposed  design 
concept  ensures  compact  representation  of  circuits  compared  with  other  formats  and  methods.  The 
following  design  algorithm  guarantees  compact  circuit  representation: 

Function  (Circuit)  <->  BDD  Model  <->  LWDD  Model  <->  Realization. 


The  LWDD  is  embedded  in  3D  X -hypercells  that  represent  circuits  in  a  3D  space.  The 
polynomial  representation  of  logical  functions  ensure  the  description  of  multi-output  functions  in  a  word- 
level  format.  The  expression  of  a  Boolean  function  f  of  n  variables  (xi,  x2,  ...  x„_i,  xn)  is 


n 

LP  =  a0  +  axxx  +  a2x2  + ...  +  an_xxn_x  +  anxn  =  a0  +  ^  a -Xj 

7=1 


.  To  perform  a  design  in  3D,  we  established  the 


following  mapping  LWDD(ao,aha2,. .  .,an_i,an)  <->  LP.  The  nodes  of  LP  correspond  to  a  Davio  expansion. 
The  LWDD  is  used  to  represent  any  m-level  circuit  with  levels  Lif  i=l,  ...,  m  with  elements  of  the 
standard  primitive  library.  Two  data  structures  defined  in  the  algebraic  form  by  a  set  of  LPs  as 
Lx  :  inputs  x.;  outputs  yxk 
L2 :  inputs  yXk ;  outputs  y2l 


L  =  \ 


Lm- 1  :  inputs  ym_2y,  outputs  ym  u. 
Lm  :  inputs  ym_hw;  outputs  ym  „ 


that  corresponds  to  LPi  =  +  fia)xJ, LPm  =  a"  +Yja"jym_XJ  , 

7=1  7=1 

or  in  the  graphic  form  by  a  set  of  LWDDs  as  LWDL\{a\,...,aln)  <->  LPX,  ...,LWDDm(al,...,ann  )<r^LPm. 


The  use  of  LWDD  is  a  significant  departure  from  existing  logic  design  tools  but  our  technology  is 
compatible  with  the  existing  software,  algorithms  and  circuit  representation  formats.  Circuit 
transformation,  format  transformation,  modular  architecture,  library  functions  over  primitives,  and  other 
features  can  be  performed.  From  the  circuit  transformation  viewpoint,  all  combinational  circuits  can  be 
represented  by  LWDD.  The  format  transformation  can  be  performed  for  circuits  defined  in  Electronic 
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Data  Interchange  Format  -  EDIF,  Berkeley  Logic  Interchange  Format  -  BLIF,  Verilog  and  other  formats. 
The  library  functions  may  have  a  library  of  LWDDs  for  multi-input  circuits,  as  well  as  libraries  of 
molecular  electronic  devices  and  Mgates.  The  important  feature  is  that  these  primitives  are  realized 
(through  logic  design),  and  should  be  fabricated  as  X -hypercells  modules.  The  reported  LWDD  provides 
modeling,  analysis,  verification,  evaluation  and  other  features. 

Arithmetic  expressions  underlying  the  design  of  LWDD  are  canonical  representations  of  logic 
functions.  They  are  alternatives  of  the  sum-of-product,  product-of-sum  and  Reed-Muller  forms  of 
representation  of  Boolean  functions.  In  LWDDs,  linear  data  structures  are  manipulated  using  the  masking 
operators.  In  particular,  LWDDs  are  obtained  by  mapping  LPs,  where  the  nodes  correspond  to  the  Davio 
expansion  and  functionalizing  vertices  to  the  coefficients  of  the  LPs.  The  LWDD  design  flow  is 

Primitive  <->  LP  Model  <->  LWDD  Model  <->  Realization. 

Any  m-level  logic  circuits  with  fixed  order  of  elements  are  uniquely  represented  by  a  system  of  m 
LWDDs  as:  Function  (Circuit)  <->  LP  Model  <->  LWDD  Model  <->  Realization.  The  proposed  concept  is 
verified  by  representing  Boolean  functions  by  hypercells  in  3D  ICs  design.  CAD  tools  for  logic  design 
must  be  based  on  the  principles  of  3D  realization  of  logic  functions  with  a  library  of  primitives.  LWDDs 
are  extended  by  embedding  the  total  solution  tree  into  the  X -hypercell  structures.  For  two  graphs  G=(V,E) 
and  H=(W,F),  we  embed  the  graph  G  into  the  graph  H.  The  information  in  the  resulting  X -hypercells  is 
subdivided  according  to  the  new  structural  properties  of  the  cell  and  the  type  of  the  embedded  tree.  The 
embedding  of  a  guest  graph  G  into  a  host  graph  His  a  one-to-one  mapping  MGV\  V(G)—>V(H),  along  with 
the  mapping  M  that  maps  an  edge  (u;v)eE(G)  to  a  path  between  MGV(u )  and  MGV(v)  in  H.  Thus,  the 
embedding  of  G  into  H  is  a  one-to-one  mapping  of  the  nodes  in  G  to  the  nodes  in  H. 

3.  1 .  3.  Information-Theoretic  Analysis 

Consider  the  input-output  information  measures  for  3D  ICs  as  shown  in  Figure  4.  One  examines 
the  information  carried  by  the  variable  and  function  f  denoted  as  /(xz)  and  1(f),  as  well  as  entropies  of 
the  input  variable  xt  and  the  output  function/  denoted  as  H(x^)  and  H(f).  The  condition  entropies  H(f\xi) 
and  H(xi\f),  relative  information  I(J\xL)  and  /(xzj/),  mutual  information  I(J,xf  as  well  as  joint  entropy  H(f,x 
are  of  our  interest. 


Input  Xi 


Figure  4.  Input-output  representation  of  MICs 

In  design,  decision  diagrams  and  decision  trees  are  used.  Information  measures  can  be  estimated 
and  evaluated.  Decision  trees  are  designed  using  the  Shannon  and  Davio  expansions  and  information 
theoretical  approaches  [12].  There  is  a  need  to  find  the  best  variable  and  expansion  for  any  node  of  the 
decision  tree  in  terms  of  information  measures  in  order  to  minimize  or  maximize  information  estimates. 
The  entropy-based  optimization  algorithm  should  generate  the  optimal  paths  in  a  decision  tree  with 
respect  to  the  design  criterion.  The  decision  tree  is  designed  by  arbitrarily  choosing  any  variable  using 
Shannon  (5),  positive  Davio  (pD )  or  negative  Davio  (nD)  expansions  for  each  node.  We  use  the 
conditional  entropy  H(J\Tree)  and  mutual  information  I(f,Tree )  as  the  performance  estimates.  The  initial 


state  of  the  design  process  is  characterized  by  the  maximum  value  for  the  conditional  entropy  H(fTree)= 
H(f).  Nodes  are  recursively  attached  to  the  decision  tree  by  minimizing  the  entropy  H(J\Tree)  and 
maximizing  the  mutual  information  I(f,Tree).  The  final  state  of  the  decision  tree  is  characterized  by 
H(J\Tree)=0  and  I(fJree)=H(f). 

The  decision  tree  design  process  is  a  recursive  decomposition  of  a  switching  function.  This 
recursive  decomposition  corresponds  to  the  expansion  of  switching  function /  with  respect  to  the  variable 
x.  The  variable  x  carries  information  that  influences  f  The  initial  state  of  the  expansion  ere  {SpDpD}  can 
be  characterized  by  the  entropy  //(/),  and  the  final  state  by  the  conditional  entropy.  In  particular,  the  cr- 
expansion  of  the  function  /  with  respect  to  the  variable  x  is  Ia(f,x)  =  H(f)  -  Ha(J\x).  The  information 
estimates  for  S ,  pD  and  nD  expansions  are  found.  If  the  design  of  a  decision  tree  is  based  on  the  S 
expansion  and  the  sum-of-products  mapping,  then,  the  information  measure  for  the  function / with  respect 
to  x  is  Hs(f\x)  =  p^x_0  -H(f\x_0)  +  p^x_x  The  information  measure  of  the  S  expansion  is  the  conditional 

entropy,  e.g.,  Hs(j\x)=H(j\x).  In  contrast,  for  the  pD  and  nD  expansions,  we  have 

H^(f\x)  =  p^.H(flx=0)+plx=l-H(flx=0®flx=1)  and  H"D(f\x)  =  pH H(flx=0  ®/M).  The 

information-theoretic  estimates  of  the  pD  and  nD  nodes,  as  compared  with  the  Shannon  expansion,  for  a 
decision  tree  design  are  MpD  =  p^\H{fv=x)-H{fv=a® f]x=x)\  and  Our 

results  justify  the  use  of  the  Davio  expansions  in  the  design  as  illustrated  for  the  benchmarking  circuits.  It 
is  found  that  the  information  estimates  are  one  of  the  criteria  to  optimize  3D  design  as  well  as  ICs.  The 
following  example  illustrates  the  concept. 

Example.  Consider  the  benchmark  Cl 7  circuit  implemented  using  3D  molecular  NAND  gates 
(mNAND)  as  reported  in  Figure  5. 


For  the  Shannon  expansion  used  in  the  design  of  decision  diagram,  from  /?(/i=0)=14/32  and 
p(/i=l)=18/32,  we  have  H% )=-(  1 4/32)log2(  1 4/32)-(  1 8/32)log2(  1 8/32)=0.99  and  Hs(f\xt)=- 
(10/32)log2(10/16H6/32)log2(6/16H4/32)log2(4/16)-(12/32)log2(12/16)=0.88.  One  also  obtains 
Hs(f 1|x2)=0.67,  Hs(f i|x3)=0.98,  H% |x4)=0.98  and  Hs(fx\x5)=0.99.  The  entropies  Hs(f2)  and  Hs(f2\x)  are 
straightforwardly  calculated.  The  output  variables  and  conditional  entropies  are  found  with  respect  to  all 
variables  if  the  positive  and  negative  Davio  expansions  are  applied.  We  have  Hs(f\ |x,-)min=0.67,  EFD(f\  Ix,-)^ 
=0.88  and TfVikOmm  =0.75. 


We  conclude  that  the  information-centered  optimization  of  ICs  design  is  formulated  as  the 
optimization  algorithm  with  respect  to  the  information  criterion  (estimates)  in  order  to  design  optimal 
decision  diagrams.  A  path  in  the  decision  tree  starts  from  a  node  and  finishes  in  a  terminal  node.  Each 
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path  corresponds  to  a  term  in  the  final  expression  for  f  The  criterion  for  choosing  the  decomposition 
variable  x  and  the  expansion  ere  {S,pD,nD}  is  that  the  conditional  entropy  of  the  function  with  respect  to 
this  variable  has  to  be  minimum,  e.g.,  min Ha(j\x). 


3.  1.  4 .  Design  of  Benchmarking  3D  ICs  and  Proof-of-Concept  Arithmetic  Logic  Unit 

The  proposed  design  of  3D  ICs  is  illustrated  and  verified  using  benchmarking  ICs.  These  ICs  are 
used  to  evaluate  and  examine  different  design  concepts.  The  design  technology  is  concurrently  supported 
by  the  representative  CAD  tools  as  reported  in  section  3.3.  A  series  of  numerical  studies  were  conducted 
to  accomplish  the  software-supported  logic  design  for  proof-of-concept  ICs  in  3D.  The  size  of  LWDDs 
generated  by  the  developed  software  is  compared  with  the  best  results  obtained  by  using  the  Decision 
Diagram  Packages,  developed  by  University  of  Colorado  and  Darmstadt  University  of  Technology,  for 
2D  VLSI  design.  The  method  was  reported  and  software  algorithms  were  tested  and  validated.  The 
number  of  nodes,  number  of  levels  and  CPU  time  (in  seconds)  required  to  design  decision  diagrams  for 
3D  ICs  are  reported  in  Table  1.  The  logic  designs  for  a  9-bit  ALU  (c5315)  and  multiplier  (c6288),  for 
which  classical  decision  diagram  may  not  be  used,  were  successfully  performed. 

Table  1.  Size  (number  of  nodes  N,  levels  L  and  CPU  time)  in  BDD,  WDD  and  LWDD 


Circuit  I/O 

BDD  WDD 

#N 

LWDD  (EDIF) 

#L  #N  CPU  Time 

LWDD  (ISCAS’85) 

#L  #N  CPU  Time 

C432  36/7 

1064  1209 

21  455  <0.01 

17  336  <0.01 

C880  60/26 

4053  4048 

35  604  <0.01 

24  605  <0.1 

c2670  233/140 

1850  3939 

40  1374  <0.1 

32  2026  <0.1 

ALU 

c5315  178/123 

1719  2504 

68  3192  <0.3 

49  4156  <0.4 

Multiplier 
c6288  32/32 

124  4320  <0.2 

124  4318  <0.5 

Our  goal  is  to  evaluate  3D  architecture  characteristics  such  as  volumistic  size,  topological 
parameters  and  performance.  We  assume:  (1)  feedforward  neural  networked  topology  with  no  feedback; 
(2)  threshold  Mgates  are  used  as  the  processing  units;  (3)  aggregated  networked  3D  X  -hypercells 
comprised  from  Mgates;  (4)  multilevel  combinational  circuits  are  used  over  the  library  of  NAND,  NOR 
and  EXOR  Mgates.  Experiments  were  conducted  for  circuits  reported  in  Tables  1  and  2.  The  space  size  is 
given  by  A,  Y  and  Z  that  result  in  the  volumistic  (V)  quantity  as  V=XxYxZ.  The  topological  characteristics 
are  evaluated  using  the  total  number  of  terminal  (#  NT)  and  intermediate  (#  A/)  nodes.  A  3D  9-bit  ALU 
(c5315)  with  178  inputs  and  123  outputs  is  implemented  using  1413  Mgates,  while  a  c6288  multiplier  (32 
inputs  and  32  outputs)  has  2327  Mgates.  As  reported  in  Table  2,  Mgates  (#G)  were  networked  and 
aggregated  (grouped)  in  3D.  The  number  of  incompletely  specified  N -hypercells  was  minimized.  The  N- 
hypercells  in  the  ith  layer  were  connected  to  the  corresponding  hypercells  in  (z-l)th  and  (z+l)th  layers. 
The  obtained  3D  architecture  has  XxYxZ  hypercells.  The  number  of  terminal  nodes  and  intermediate 
nodes  are  3750  and  2813  for  a  9-bit  ALU,  while,  for  a  multiplier,  we  have  9248  and  6916  nodes.  To 
combine  all  layers,  more  than  10,000  connections  were  generated.  The  design  in  3D  was  performed 
within  0.4  seconds.  Other  metrics  to  evaluate  the  designed  3D  ICs  were  used.  That  is,  in  addition  to 
conventional  parameters  (diameter,  dilation  cost,  expansion,  load,  etc.),  we  used  the  number  of  variables 
in  the  logic  function  described  by  N -hypercells,  number  of  links,  fan-out  of  the  intermediate  nodes, 
statistics,  etc. 
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Table  2.  Experimental  results  for  3D  ICs 


Circuit 

I/O 

#G 

Space  Size 
#X  #Y 

#z 

Nodes  and  Connections 
#Nr  #  /V/  CPU  Time 

c432 

36/7 

126 

66 

64 

66 

2022 

1896 

<0.01 

c880 

60/26 

130 

70 

72 

70 

612 

482 

<0.1 

c2670 

233/140 

828 

82 

80 

78 

3594 

2766 

<0.1 

ALU 

c5315 

178/123 

1413 

138 

132 

126 

3750 

2813 

<0.4 

Multiplier 
c6288  32/32 

2327 

248 

248 

244 

9246 

6916 

<0.5 

Comparison  of  2D  and  3D  designs  for  a  9-bit  ALU.  The  high-level  model  is  documented  in 
Figure  6.  The  studied  ALU  performs  arithmetic  and  logic  operations  simultaneously  on  two  9-bit  input 
data  words,  as  well  as  computes  the  parity  of  the  results.  Conventional  2D  logic  design  of  a  9-bit  ALU 
(c5315)  with  178  inputs  and  123  outputs  results  in  2406  gates.  In  contrast,  the  proposed  design  (as  an 
envisioned  SLSI)  leads  to  1413  Mgates  that  are  networked  and  aggregated  in  3D. 
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3.  2.  X -HYPERCELL  DESIGN 


The  binary  tree  is  a  networked  architecture  that  carries  information  about  dual  connections  of 
each  node.  The  binary  tree  also  carries  information  about  functionality  of  the  logic  circuit  and  topology. 
The  nodes  of  the  binary  tree  are  associated  with  the  Shannon  and  Davio  expansions  with  respect  to  each 
variable  and  coordinate  in  3D.  A  node  in  the  binary  decision  tree  realizes  the  Shannon  decomposition 
/  =  x./0  ©  x.J\ ,  where  fQ  =  /|  q  and  fl=  f\  _1  for  all  variables  in  f  Thus,  each  node  realizes  the 

Shannon  expansion,  and  the  nodes  are  distributed  over  levels.  The  classical  hypercube  contains  2n  nodes, 

n 

while  the  X -hypercells  has  2n  +  ^2 n~lC™  nodes  in  order  to  design  technology-dependent  biomolecular 

z=0 

ICs  in  3D.  The  X -hypercell  consists  of  terminal  nodes,  intermediate  nodes  and  roots.  The  most 
straightforward  X -hypercell  implementation  is  provided  by  using  molecular  multiplexers  (mMUX). 
However,  other  Mgates  can  be  utilized  as  well.  The  design  steps  are: 

Step  1 :  Connect  the  terminal  node  with  the  intermediate  nodes; 

Step  2:  Connect  the  root  with  two  intermediate  nodes  located  symmetrically  on  the  opposite  faces; 
Step  3:  Pattern  the  terminal  and  intermediate  nodes  on  the  opposite  faces  and  connect  them  via 
the  root. 

The  nanobiocomputing  architectures  can  be  designed  and  implemented  as  aggregated  X- 
hypercells  engineered  (fabricated)  from  organic  or  bio  molecules.  Figure  7. a  reports  a  3D  X -hypercell 
which  comprises  from  a  two-to-one  molecular  multiplexer  (mMUX).  A  Molecular  multiplexer  and  other 
Mgates  (Minverter,  MAND,  MNAND,  MOR,  MNOR  and  MXOR),  designed  by  the  Microsystems  and 
Nanotechnologies  team  using  cyclic  carbon-based  molecules,  are  reported  in  Figure  7.b. 

There  are  several  methods  for  representing  logic  functions.  A  X -hypercell,  proposed  as  a  core 
solution,  is  a  homogeneous  aggregated  assembly  for  massive  super-high-performance  parallel  computing. 
We  apply  the  enhanced  switching  theory  integrated  with  a  novel  logic  design  concept  [12].  In  the  design, 
the  graph-based  data  structures  and  3D  circuit  topology  are  utilized.  The  X -hypercell  is  a  topological 
representation  of  a  switching  function  by  a  ^-dimensional  graph.  In  particular,  the  switching  function  /  is 
given  as 

Coefficient 

2n _ 1  II 

1  /  \  Form  of  Switching  Function 

/  =>  L  K;  (*?...*?)=>  fF 

i= 0 

ft 

Operation 

The  data  structure  is  described  in  matrix  form  using  the  truth  vector  F  of  a  given  switching 
function / as  well  as  the  vector  of  coefficients  K.  The  logic  operations  are  represented  by  L.  X -hypercells 
compute  /  and  complex  switching  functions  are  implemented  by  hypercells.  For  example,  Figure  8 
reports  a  X -hypercell  to  implement  /  =  x1x2  v  x^x2  v  x1x2x3 .  From  the  technology-specific  viewpoints,  we 

propose  a  X-hypercell  solution  that  employs  Mgates  as  documented  in  Figures  7  and  8.  This  3D  solution 
coherently  maps  the  device-  and  gate-level  consideration  by  aggregating  and  networking  Mgates.  By 
utilizing  the  root  and  intermediate  nodes  at  the  edges,  as  shown  in  Figure  8,  the  reported  X -hypercells 
implement  switching  functions  / of  arbitrary  complexity. 
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Figure  7.  a)  K -hypercell  utilizing  two-to-one  multiplexer  (mMUX); 
b)  Mgates  (Minverter,  MAND,  MNAND,  MOR,  MNOR,  MXOR  and  mMUX)  which  can  be 
utilized  in  the  design  of  K -hypercells 


The  logic  design  in  spatial  dimensions  was  based  on  the  advanced  methods  and  enhanced 
data  structures  to  satisfy  a  3D  topology.  The  appropriate  data  structure  of  logic  functions  and  methods  of 
embedding  this  structure  into  X -hypercells  were  found.  The  three-step-solution  in  a  logic  functions’ 
manipulation  in  order  to  change  the  carrier  of  information  from  the  algebraic  form  (logic  equation)  to  the 
hypercell  structure  is: 

Step  1\  Logic  function  is  transformed  to  the  appropriate  algebraic  form  (Reed-Muller,  arithmetic 
or  word-level  in  matrix  or  algebraic  representation); 

Step  2\  Derived  algebraic  form  is  converted  to  the  graphical  form  (decision  tree,  decision 
diagram  or  logic  network); 

Step  3:  Obtained  graphical  form  is  embedded,  and  technologically  implemented,  by  X- 
hypercells. 
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Figure  8.  K -hypercell  that  implement  function / 
The  design  is  expressed  as 

Logic  Function  <=>  Graph  <=>  N  -  Hypercell  Structure  . 

Step  1  Step  2  Step  3 


The  proposed  procedure  results  in: 

•  Algebraic  representation  and  robust  rules  of  manipulation  for  complex  switching  logic 
functions; 

•  Matrix  representation  and  sound  manipulation; 

•  Consistency  of  logic  relationships  for  variables  and  functions  from  the  spectral  theory 
viewpoint  due  to  the  use  of  matrix  algebra; 

•  Graph-based  representation  is  found  using  decision  trees,  decision  diagrams  and  logical 
networks; 

•  Data  structures  are  embedded  into  X -hypercells. 


3.  3.  THREE-DIMENSIONAL  NANOBIOCOMPUTING  CAD 


Fully  featured  technology-dependent  CAD  for  SLSI  design  is  a  significant  task  that  is  not  within 
the  scope  of  this  contract.  The  representative  CAD  tools  and  software  solutions  were  developed  in  order 
to  demonstrate  the  3D  design  feasibility.  In  particular,  we  developed  a  Linear  Decision  Diagram  (LDD) 
Package  and  Windows-Based  3D  ICs  Interactive  Toolbox  (WindowIT)  using  C  and  OpenGL.  We  plan  to 
integrate  both  packages.  The  compatibility  with  hardware  description  languages  is  important.  To  be 
compatible  with  conventional  logical  programming  tools,  three  netlist  formats  (Electronic  Data 
Interchange  Format  -  EDIF,  International  Symposium  on  Circuits  and  Systems  -  ISCAS,  and  Berkeley 
Logic  Interchange  -  BLIF)  are  used  and  embedded  in  the  LDD  Package.  The  methods  reported  in  sections 
3.1  and  3.2  were  used  in  the  software  development.  A  coherent  data-intensive  analysis  of  SLSI-designed 
ICs  can  be  performed  by  a  completely-featured  CAD  as  molecular  bioelectronics  technology  is 
developed.  We  can  examine  baseline  characteristics  and  prove  the  efficiency,  robustness,  feasibility  and 
other  significant  advantages  of  the  proposed  technology. 
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The  representative  components  of  the  developed  CAD  prototype  (LDD  Package  and  WindowIT), 
allow  one  to  carry-out  the  design  of  ICs  in  3D,  as  well  as  to  examine  major  performance  estimates.  The 
software,  developed  by  Microsystems  and  Nanotechnologies ,  was  successfully  tested  and  verified.  The 
LDD  Package,  based  on  our  synthesis  concept  as  reported  in  sections  3.1  and  3.2,  features: 

-  new  sound  design  technology  and  methods  for  3D  ICs; 

-  synthesis  and  partitioning  linear  decision  diagrams  for  given  functions  or  circuits; 

-  spectral  representation  of  logic  functions; 

-  circuit  verification; 

-  compact  format  for  rapid-prototyping; 

-  compressed  representation  of  complex  logic  networks. 

Results  in  design  and  visualization  of  3D  ICs  are  reported  in  Figure  9  displaying  the  Command 
Windows  data  for  the  LDD  Package.  It  should  be  emphasized  that  the  results  of  the  3D  design  of  a  cl7 
circuit  as  well  as  a  9-bit  ALU  (c5315)  are  displayed.  The  WindowIT  environment  is  documented  in  Figure 
10.  As  were  reported  in  sections  3.1  and  3.2,  we  were  able  to  successfully  carry-out  and  verify  3D  design 
for  benchmarking  ICs. 
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LDDPackage 


I  Auio  21  nNffll  ES]  iff[g  a| 
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build_ldd_l  <bll> 
help 
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print_inputs  <pi) 
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print_ldd_stats  Cpls> 
print_network  CpiD 
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print_outputs  <po> 
quit 

read_blif  filename  (rb> 
read_edif  filename  (re> 
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read_ldd_a  filename  Crla> 
read_ldd_l  filename  (rll) 
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simulate_ldd  <sl> 
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Figure  9.  Design  of  3D  ICs  using  LDD  Package  that  performs  integrated  design,  aggregation,  analysis, 

optimization  and  evaluation 
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Figure  10.  Design  of  3D  ICs  using  interactive  WindowIT  that  performs  integration,  aggregation 

and  visualization 
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3.  4.  DNA  DERIVATIVE  TRANSISTOR  FOR  X -HYPERCELLS 


Sequence-dependent  self-assembled  DNA  and  templated  protein  synthesis  can  be  used  to  build 
patterned  two-  and  three-dimensional  structures  with  the  desired  geometrical  topologies  [13,  14].  Some 
inorganic,  organic,  hybrid  and  biomolecules  which  exhibit  desired  electronic  properties  can  be  utilized  in 
a  new  generation  of  electronic  nanodevices  and  in  nano-bio-electronic  architectures  [2,  4,  15-17]. 
Complex  3D  information  processing,  computing  and  memory  architectures  are  envisioned  to  be  designed 
and  fabricated  using  self-assembled  biomolecules  (including  DNA  and  proteins)  utilized  as  enhanced- 
functionality  multi-terminal  interconnected  electronic  nanobiodevices.  It  is  important  to  design  and 
analyze  functional  high-performance  bioelectronic  devices  while  taking  basic  phenomena  and  effects  in 
biomolecule-interconnect  complexes  into  account.  Insulating,  conducting  and  semiconducting  properties 
have  been  observed  in  the  contact-biomolecule-contact  and  inter-biomolecular  complexes  [18-20].  DNA 
synthesis,  sequence,  length,  environment,  alignment  and  other  factors  influence  the  electronic 
characteristics.  In  DNA,  various  experiments  to  examine  the  charge  and  electron  transport,  I-V 
characteristics,  conductance  and  other  electronic  properties  were  performed  in  [18-20]  with  limited 
success.  Figure  11. a  illustrates  DNA  attached  to  two  electrodes  [19].  The  electrode  (source  and  drain) 
electrochemical  potentials  are  denoted  as  VFs  and  VFd,  while  Es  and  Ed  are  the  self-energy  functions  of  the 
source  and  drain.  Electrochemical  potentials  VFs  and  VFd  vary.  There  is  no  electron  transport  if  the  system 
is  in  equilibrium,  e.g.,  VFs  =  VFd.  The  highest  occupied  molecular  orbitals  (HOMO)  and  lowest 
unoccupied  molecular  orbitals,  (LUMO),  as  well  as  the  Fermi  level,  are  illustrated.  Depending  on  the 
HOMO  and  LUMO  levels,  as  well  as  £>,  the  electron  transport  takes  place  through  the  particular  orbitals. 
Using  broadening  energies  EBs  and  EBd ,  the  electron  transfer  rates  are  EBs/h  and  EBd/h.  Guanine  quartets 
(tubular  sequences  of  G  tetramers)  were  reported  [19].  Tetramers  (G4)  are  the  building  blocks  of  a 
quadruple-helix  forming  G4-DNA.  The  hydrogen-bonded  guanines  arrange  in  a  2.3  nm  diameter 
configuration.  The  G4-DNA-based  multi-terminal  electronic  nanodevice  is  illustrated  in  Figure  ll.b. 
Though  G4-DNA  aggregates  were  fabricated,  the  electronic  characteristics  have  not  been  experimentally 
examined  due  to  unsolved  interconnect,  manipulation  and  characterization  difficulties. 


(a)  (b) 

Figure  1 1.  a)  Double-stranded  DNA  with  two  electrodes; 
b)  G4-DNA-based  multi-terminal  electronic  nanodevice  and  3D  nanobioICs  topology 


Due  to  the  limited  experimental  data  for  the  electron  transport  and  I-V  characteristics  in 
biomolecules,  we  examined  a  hybrid  CMOS -technology  device  for  which  the  experimental  data  is 
available  [21,  22].  Prototypes  of  hybrid  bioelectronic  devices  were  built  by  utilizing  enhanced  CMOS 
technology.  The  electronic  behavior  and  I-V  characteristics  in  the  planar  topology  DNA-based  transistor 
(DNAt),  as  reported  in  Figure  12. a,  were  studied  merging  experimental  and  theoretical  results.  The 
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studied  DNAT  is  equivalent,  to  some  extent,  to  the  silicon  technology  CMOS  field-effect  transistors 
(FET).  For  the  planar  DNAT,  the  experimental  I-V  characteristics  were  measured  as  documented  in 
Figure  12.b. 


Experimental  I-  V  Characteristics  of  DNA-Based  FET 
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Figure  12.  DNAT  and  experimental  I-V  characteristics 


The  fabrication  of  this  DNAT  is  performed  utilizing  CMOS  technology  and  depositing  A-DNA  on 
silicon  oxide  with  the  adsorption-based  source-DNA-drain  interconnect.  It  is  reported  that  an  A-DNA 
bundle  of  poly(dA)-poly(dT)  is  a  ft-type  semiconductor,  while  poly(dG)-poly(dC)  is  a  /7-type  [22].  DNA 
derivatives  may  ensure  self-assembled  and  ordered  engineered  bundles,  layers  and  ropes  which  are 
organized  and  interconnected  due  to  specific  reactivity  of  molecules’  functional  groups  with  affinity  to 
distinct  surfaces  and  molecules.  The  guanine  (DNA  nitrogenous  base)  exhibits  a  low  oxidation  potential 
that  results  in  electron  transport  and  self-assembly  in  liquid  and  solid  media  [23].  DNAT  with  self- 
organized  H-bonded  DNA  bases  were  studied  in  [21-23]  utilizing  a  polymerized  derivative  of 
deoxyguanosine  poly(dG)„.  A  single  modified  DNA  base,  e.g.,  poly(dG)„,  is  deposited  on  the  silicon 
oxide  between  source  and  drain  as  reported  in  Figure  12. a  forming  the  channels.  The  experimental  I-V 
characteristics  are  documented  for  a  different  voltage  uGSr,  which  controls  the  electron  flow  in  the  source- 
drain  channel. 


The  modeling  is  performed  using  the  Green’s  function  formalism  [4,  24-27].  The  Hamiltonian 
(H)  for  the  open  contact  (left  and  right  electron  reservoirs)  -  biomolecule  (bridge)  complex  is  constructed 
by  using  the  biomolecule  HM,  tunneling  HT,  contact  (junction)  Hc  and  external  HE  Hamiltonians.  That  is, 
we  have  H  =  HM  +  HT  +  Hc  +  HE.  To  derive  the  I-V  characteristics  of  the  biomolecular  complex 
(electrode-biomolecule-electrode),  we  apply  Keldysh’s  nonequilibrium  Green  function  concept.  Green’s 
function  is  a  wave  function  of  energies  at  r  resulting  from  an  excitation  applied  at  rE.  We  applied  the 
retarded  Green’s  function  G  that  represents  the  response  of  the  biomolecule-contact  complex.  Hence, 
(E-H)G(r,rE)  =  S( r-rE).  The  boundary  conditions  must  be  satisfied  for  the  transport  and  Poisson 
equations.  The  molecular  Hamiltonian  of  the  isolated  system  and  the  complex  self-energy  functions  are 
used  to  examine  the  finite  biomolecule  complex,  instead  of  the  single  energy  potential  and  broadening 
energies.  In  the  matrix  notations,  using  the  overlap  identity  matrix  S,  the  Green’s  function  is  found  to 
be [G’( E)]  =  (zTS’]  —  [//]— [Vsc ]  -  Z  [£;  ])'  >  where  [£’(.]  =  [5,.][GI.][5'*];  S,  is  the  geometry-dependent 

coupling  (contact)  matrix  between  the  surface  of  the  source/drain  and  biomolecule  (for  source  and  drain 
Ss  and  Sd  are  the  geometry-dependent  coupling  matrices  between  the  surface  of  the  source  and  drain  and 
the  device).  The  imaginary  non-Hermitian  self-energy  functions  of  the  source  and  drain  (first  and  second 
electron  reservoirs)  Es  and  Ed  are  [Es]  =  [*SJ[GJ[*S'*]  and  [Ed]  =  [Sd][Gd][S*d],  where  Gs  and  Gd  are  the 
surface  Green’s  functions  for  the  source  and  drain  that  are  found  applying  the  recursive  methods. 
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Taking  note  of  the  Green’s  function,  the  density  of  state  D{E)  is  found  using  G{E)  as 
D(E )  =  ~^\m{G(E)} .  The  spectral  function  A(E)  is  the  anti-Hermitian  term  of  the  Green’s  function,  and 

A{E )  =  i[G(E )  -  G* (£)]  =  -2  Im[G(E)] .  One  obtains  D(E)  =  E  tr[^(£,)iS] ,  where  tr  is  the  trace 

operator,  and  S  is  the  identity  overlap  matrix  for  orthogonal  basis  functions.  Using  the  broadening  energy 
functions  EBs  and  EBd ,  one  obtains  the  source  and  drain  spectral  functions 
[Am  =  [G(E)][EBs(E)][G\E)]  ,  [A.m  =  [G(E)][EBd(E)][G\E)] . 


For  multi-terminal  electronic  nanobiodevices,  the  nonequilibrium  charge  density  matrix  is 

-.00  ..00 

[pm  =  —  \  Z  f(Evk,VFij)[AjmdE  =  —  J  Z  f(El/km,)[G(E)][EBii(E)][G'(E)]dE, 

_oo  k,ies,jed  _oo  k,ies,jed 

where  VFij-  are  the  potentials,  e.g.,  VFs  and  VFd  are  the  source  and  drain  potentials  related  to  the  Fermi 
levels;  f[Ev,VFiJ)  are  the  Fermi-Dirac  functions,  for  example,  J[Ev,VFs)  and  f[Ev,VF(j)\  Ev  is  the  single 
energy  potential  that  depends  on  the  charge  density  p(E )  ,  Ev  =  Ev0  +  p(E)  .  Utilizing  the  transmission 

matrix  T(E)  =  tr [EBsG(E)EBdG*(E)]  and  taking  note  of  the  Green’s  function  and  broadening,  the  current  is 

+00 

found  to  be  /t=-U tr[EBsG(E)EBdG* (E)]  Z REn,VFiJ)dE  • 

— oo  k,iGS,j&d 


The  contact-biomolecule-contact  complex  attains  equilibrium  with  the  Fermi  level.  The  current- 
voltage  characteristics  are  calculated.  The  solutions  and  I-V  characteristics  are  found  carrying  out  robust 
heterogeneous  MATLAB  simulations  with  data-intensive  analysis  features.  The  non-equilibrium  charge 

oo 

density  is  p(E)  =  —  f  [f(Ev,VFs)G(E)EBsG\E)  +  f(Ev,VFd)G(E)EBdG\E)\JE  ,  where  J[Ev,VFs)  and 
2n  J 


f(Ev,VFd)  are  the  Fermi-Dirac  functions,  f(Ev, VFs)-- 


e,-vFs  ■> 

-1 

(  Ey~VFd  \ 

1  +  e  kT 

’  f(Ev,VFd)  = 

l  +  e  kT 

J 

v  J 

.  Hence,  the 


ry  +00 

current  between  two  junctions  (contacts)  is  /  =  —  J  tr  [EBsG(E)EBdG*  (E)\f(Ev  ,VFs)-  f{Ev ,  VFd  )}lE  . 

^  -00 

For  the  numerical  studies,  the  LUMO  and  HOMO  energy  gaps,  broadening  energies,  Fermi 
energy  £>,  as  well  as  other  data  are  used.  It  should  be  emphasized  that  the  polymerized  derivative  of 
deoxyguanosine  poly(dG)„  forms  non-uniform  and  non-perfect  channels  of  different  length  with  distinct 
binding  to  the  source  and  drain.  Furthermore,  the  electrostatic  field  that  affects  the  electron  transport  is 
non-uniform.  Correspondingly,  in  general,  one  may  not  be  able  to  achieve  a  precise  matching  of  the 
experimental  and  analytical  (numerical)  results.  The  energy  levels  of  the  molecular  orbitals,  examined  for 
the  free  biomolecule,  thiol-biomolecule,  as  well  as  the  Au-S-biomolecule-S-Au  complex,  are  different. 
In  analysis  of  electronic  characteristics,  one  considers  the  biomolecule  orbitals;  3s  sulfur  orbitals,  as  well 
as  5d,  6s  and  6p  gold  orbitals.  The  HOMO  and  LUMO  of  the  Au  cluster  shift  towards  lower  energies 
upon  interaction  with  the  functionalized  biomolecule.  The  contribution  of  different  orbitals  to  the  electron 
transport  varies.  Two  occupied  orbitals  result  in  a  strong  coupling  (hybridization)  between  gold,  sulfur 
and  the  ^--conjugated  molecule.  These  orbitals  are  of  great  importance  for  electron  transport  between 
contacts  and  the  functionalized  molecule.  The  HOMO-LUMO  gap  is  1.75  eV.  However,  other 
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LUMOs  and  HOMOs  contribute  to  the  electron  transport.  Simulation  results  are  reported  in  Figures  13  to 
15  (dashed  lines)  for  different  coupling  and  Fermi  energies  of  the  gold  cluster,  assuming  a  single 
biomolecule  energy  level  Ev.  One  should  integrate  multiple  HOMO  and  LUMO  states.  Therefore,  in 
addition  to  single-state  modeling  (numerical  results  are  reported  in  Figures  13  to  15  by  dashed  lines), 
modeling  and  simulation  with  multi-level  discrete  energy  levels  were  performed.  The  results  are  shown  in 
Figures  13  to  15  as  a  solid  line.  The  value  of  Evo  is  Evo  =  (ELUMO-  1.5)  eV.  Figure  16  reports  the  I-V  and 
G-V  characteristics  if  Evo  is  close  to  LUMO,  and,  in  particular,  Evo  =  ( ELUMO  -  0.2)  eV.  The  results 
qualitatively  explain  and  quantitatively  match  the  experimental  measurements  of  the  I-V  characteristics 
with  distinct  coupling  energies,  Fermi  level  and  initial  single  energy  potential  Evo.  The  results  reported 
provide  evidence  that  high-fidelity  modeling  and  data-intensive  analysis  can  be  performed  using  the  non¬ 
equilibrium  Green  function  method.  Furthermore,  the  I-V  characteristics  can  be  controlled  (shaped)  and 
varied  utilizing  organic  synthesis  and  design. 
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Figure  13.  The  I-V  and  G-V  characteristics  of  DNAT  with  EBL  =  EBR  =  0.1  eV 
a)  EF  =  ( Elumo  -  0.4)  eV ;  b )  EF=  ( ELUMO  -  0.8)  eV ;  c)  EF  =  ( ELUM0  -  1 .2)  eV 
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Figure  14.  The  I-V  and  G-V  characteristics  of  DNAT  with  EBL  =  EBR  =  0.2  eV 
a)  EF  =  (Elumo  -  0.4)  eV ;  b)  EF  =  (ELUMO  -  0.8)  eV ;  c)  EF  =  ( ELUM0  -  1 .2)  eV 
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Figure  15.  The  I-V  and  G-V  characteristics  of  DNAT  with  EBL  =  0.1  eV,  EBR  =  0.3  eV 
a)  EF  =  ( ELjjMo  —  0.8)  eV ;  b)  EF  =  (ELUMO  —  1 .2)  eV 
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Figure  16.  The  I-V  and  G-V  characteristics  of  DNAT  with  EBL=  0.1  eV,  EBR=  0.3  eV,  Evo=  (ELUMO- 0.2)  eV 
a)  EF  =  ( Elumo  -  0.4)  eV;  b)  EF  =  (ELUMO  -  0.8)  eV 


The  experimental  and  numerical  results  reported  show  a  significant  variation  of  electron  transport 
in  the  DNAT  depending  on  potential  barriers,  charge  distribution,  quantization  of  energy  levels,  quantum 
interference,  binding  to  the  electrodes  and  other  effects.  The  obtained  I-V  characteristics  through 
modeling  and  simulation  match  the  experimental  data.  The  numerical  and  experimental  data  demonstrates 
the  validity  of  the  results.  Therefore,  we  achieved  accuracy  and  coherence  in  the  ability  to  analyze 
biomolecular  electronic  devices. 
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4.  Conclusions 


As  documented  in  this  report,  Microsystems  and  Nanotechnologies  has  successfully  completed  all 
Tasks  outlined  in  the  Statement  of  Work.  A  novel  design  concept  for  3D  ICs  and  computing  architectures 
was  developed,  verified  and  demonstrated.  The  proposed  design  is  supported  by  novel  software  and 
envisioned  hardware  technologies.  For  complex  ICs,  we  performed  the  design  by  utilizing  X -hypercells. 
At  the  system  level,  benchmarking  ICs  (that  have  been  widely  used  to  evaluate  2D  VLSI/ULSI  design) 
were  examined  as  proof-of-concept  3D  ICs  for  the  envisioned  biomolecular  computing  architectures.  At 
the  device  level,  we  examined  a  DNAT,  as  a  possible  proof-of-concept  three-terminal  hybrid  bioelectronic 
device.  The  aggregated  X  -hypercells  form  neuronal  3D  X -hypertopologies,  thereby  implementing  ICs 
and  nanobiocomputing  architectures.  The  representative  CAD  components  were  developed  in  order  to 
coherently  support  the  reported  SLSI  design  for  3D  ICs.  Microsystems  and  Nanotechnologies  has  proposed 
a  revolutionary  3D3  (Hardware-Software-Nanotechnol°gy)  technology  that  ensures  super-high-performance  computing 
and  information  processing.  We  focused  on  development  of  a  feasible,  practical,  affordable  and  superior 
technology  for  massive  parallel  distributed  computing  and  information  processing  for  Air  Force  systems. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


2D 

-  Two-dimensional 

3D 

-  Three-dimensional 

ALU 

-  Arithmetic  logic  unit 

BJT 

-  Bipolar  junction  transistor 

CAD 

-  Computer-aided-design 

CMOS 

-  Complimentary  metal-oxide  semiconductor 

DD 

-  Decision  diagram 

DNAt 

-  DNA  transistor 

FET 

-  Field-effect  transistor 

G-V 

-  Conductance-voltage  characteristic 

HOMO 

-  Highest  occupied  molecular  orbital 

K -hypercells 

-  Neuronal  hypercell 

ICs 

-  Integrated  circuit 

MICs 

-  Molecular  integrated  circuit 

I-V 

-  Current-voltage  characteristic 

Mgates 

-  Molecular  gate 

LDD 

-  Linear  decision  diagrams 

LP 

-  Linear  arithmetic  polynomial 

LUMO 

-  Lowest  unoccupied  molecular  orbital 

LWDD 

-  Linear  word-level  decision  diagrams 

SLSI 

-  Super-large-scale  integration 

VLSI 

-  Very-large-scale  integration 

ULSI 

-  Ultra-large-scale  integration 

29 


